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Abstract 

We consider hyperbolic equations with time-dependent coefficients and develop 
an abstract framework to derive the asymptotic behaviour of the representation 
of solutions for large times. We are dealing with generic situations where the 
large-time asymptotics is of hyperbolic type. 

Our approach is based on diagonalisation procedures combined with asymp¬ 
totic integration arguments. 


1. Introduction 

In a first part [20] of the present paper written jointly with M. Ruzhansky 
we developed tools to deal with generic hyperbolic systems with t-dependent 
coefficients based on natural assumptions and to derive energy and dispersive 
type estimates, but mainly focussed on the treatment of high frequencies and 
on estimating resulting Fourier integrals. The treatment of small frequencies 
in [20] was quite basic and with assumptions tailored towards a simple and 
rough estimate. Here we will complement and extend this paper by giving an 
alternative and more systematic approach to the treatment of small frequencies. 

In [15] a new method was developed by the author jointly with W. Nunes 
do Nascimento to describe the large-time behaviour of solutions to hyperbolic 
equations with time-dependent coefficients and to characterise it in terms of 
large-time principal symbols, see Sections 2.1.1 and 2.4 for an explanation. This 
technique will be developed further in this paper and in particular applied to 
a broad range of hyperbolic systems and equations with time-dependent coeffi¬ 
cients. 

Before giving model cases and admissible examples for our treatment, we 
will recall a bit of the history of the underlying problem and the mathematical 
background. Dispersive type estimates are an essential tool to study non-linear 
equations and to prove stability and scattering results, but they are also inter¬ 
esting on their own right due to their relation to Fourier restriction theorems 
and to decay properties of Fourier transforms of surface carried measures. For 
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the wave equation, dispersive estimates go back to the work P. Brenner [1] and 
R.S. Strichartz [23]. Higher order problems were first considered by M. Sugimoto 
[24, 25] based on refined stationary phase estimates for slowness surfaces with 
critical points. This was the starting point for a systematic treatment of decay 
properties of solutions to general scalar higher order equations with constant co¬ 
efficients and to hyperbolic systems with constant coefficients as carried out by 
M. Ruzhansky and J.G. Smith [19]. Equations of second order and with time- 
dependent coefficients have been intensively studied in the literature, examples 
include the work of M. Reissig, K. Yagdjian and J.G. Smith [16, 17, 18] or the au¬ 
thors own work on dissipative wave equations [27, 10, 28]. There are also known 
results on hyperbolic systems with time-dependent coefficients or equations of 
higher order, see e.g., the work of M. d’Abbico, S. Lucente, G. Taglialatela, 
M. Reissig and M. Ebert [5, 6, 7], and approaches based on asymptotic in¬ 
tegration techniques by T. Matsuyama and M. Ruzhansky [12, 13, 14] with 
applications to the well-posedness of Kirchhoff equations. 

All these approaches are tailored to specific situations. This has the ad¬ 
vantage of obtaining very precise conditions on coefficients and gives also the 
possibility of slight improvements of decay rates and / or assumptions on initial 
data. What we aim for in this paper and to contrast these earlier investigations 
is a general treatment of most of these cases allowing for a more structural 
understanding of the properties of such equations. 

The structure of this paper is as follows. First we will explain several model 
cases and examples which can be treated by our methods. This will be done in 
Section 2 combined with the precise description of our basic assumptions. After 
this. Section 3 will deal with the construction of representation of solutions and 
Section 4 provides the resulting energy and dispersive type estimates. Finally, 
Section 5 comes back to the examples and explains what kind of results one can 
obtain in such particular circumstances. 

In order to simplify notation, we write f ^ g for two functions / and g if 
there exists a constant C such that / < Cg uniform in the arguments. We also 
write / >< 3 if / .$ 5 and g f ■ If A: is a parameter or one of the arguments, we 
write / <7 to emphasise that the constant depends on the parameter k. We 

use a similar notation for asymptotic inequalities, f(t) < g{t), t —>■ oo, means 
that there exists a number to such that f ^ g on {t : t > fo}- 


2. Model cases and basic assumptions 


2.1. Notation and motivating examples 

We recall notation used in [20] and [22]. For A G R we denote by 




feC° 






( 2 . 1 ) 


the set of smooth functions satisfying a suitable symbol like behaviour. 
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We will treat the T-classes as classes of coefficient functions for hyperbolic 
problems. To motivate the later treatment we consider first some examples. 
They all fit into the general scheme 


DtU = Ait,D,)U, U{0,-) = Uo, (2.2) 

of systems with time-dependent matrix valued Fourier multipliers A{t,T)x) sat¬ 
isfying uniform strict hyperbolicity assumptions combined a suitable estimates 
for the t-dependence. Precise conditions are given in Section 2.4. 

2.1.1. Differential hyperbolic systems 

For coefficient matrices Aj G T{0}(8)C'^^‘^ and B G we consider 

the Cauchy problem 


DtU = Y,Aj{t)D,.U+ B(t)U, UiO,-) = Uo, (2.3) 

where D = —i9, x G M" and t > 0. In order for this system to be uniformly 
strictly hyperbolic we assume that the matrix-valued hyperbolic principal symbol 

n 

Gl(t,e) = ^A,(t)$, (2.4) 

i=i 

has for f G = {^ G K." : |CI = 1} t > 0 uniformly distinct real 

eigenvalues. We assume further the following technical-looking condition on the 
existence of a uniformly bounded and invertible diagonaliser M{t,ff) of A{t,f) 
such that 

Im f diag(M"i(r,^)B(T,0M(T,0-h (D^M"^(T,C))M(T,C))dT 
J S 

X (log(l+t)-log(l + s))I (2.5) 

holds true uniformly in t, s and f. It turns out that the leading terms of the 
integrand are independent of the particular choice of the diagonaliser M and 
should be considered as a hyperbolic subprincipal symbol of the system. It is 
shown in [22, Thm. 4.4] that the boundedness of the left-hand side is equivalent 
to a form of generalised energy conservation for high frequencies. See also [11] 
for a simpler case. The condition can be reduced to spectral conditions on the 
matrix B if the all the Aj are self-adjoint. 

Both conditions are sufficient to give a full description of the behaviour of the 
spatial Fourier transform U{t,f,) of solutions for (I -|- t)|^| > I. The large-time 
behaviour of solutions in the remaining part of the extended phase depends on 
a second principal symbol. Multiplying equation (2.3) by t gives the Fuchs-type 
equation 

n 

®tC/ = ^a,(t)D,,C/ + tB(t)C/, (2.6) 
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and we assume that 


r°° dr 

lim = Boo G / ||TB(r)-Boo|r— <oo (2.7) 

t—^oo T 

for some diagonalisable matrix B^o and a number cr > 1. The eigenvalues of 
Boo determine the large-time behaviour of solutions for (1 -|-t)|^| < 1. We refer 
to the matrix Boo (or rather to the pair consisting of the eigenvalues and the 
eigenvectors, see Section 2.4) as the large-time principal symbol of the hyperbolic 
system (2.3). Together with the constants in the estimates (2.5) it determines 
the large-time behaviour for all frequencies. 


2.1.2. Hyperbolic systems of higher order 

Hyperbolic systems of higher order arise for example when studying wave 
phenomena in elasticity. We consider only a particular model 

n n 

d2[/= + Bo{t)DtU + J2Bj{t)T)^,U + C{t)U, (2.8) 

*j=i i=i 

B(0,-) = C/o, DtB(0, •) = [/!, (2.9) 

with coefficients Aij G T{0} 0 Bj G T{1} 0 and C G T{2} 0 
subject to natural conditions. First, we assume that the matrix-valued function 

n 

^ ( 2 . 10 ) 

* j=i 


has distinct and strictly positive real eigenvalues uniform with respect to t > 0 
and f G We define the auxiliary function (later we need a smoothed 

version of this) 


h{t,£) 


ICI, 


N 


{i+tm<N. 


( 2 . 11 ) 


It allows to rewrite the second-order system (2.8) as system of first order in 
V = {BtU,h{t,B,)U)^, 


BtV=( (Gl(t,D,) + E;=ii?j(t)D,,+C'(t))h(t,D,)-i 

Vh(t,D,)I (Dth(t,D,))h(t,D,)-il 

( 2 . 12 ) 

This system of first order is pseudo-differential. Assumption (2.10) implies its 
uniform strict hyperbolicity, the hyperbolic subprincipal symbol involves A(t, f) 
in combination with the coefficients Bj{t), while the large-time principal symbol 
is determined by the matrices Bo(t) and C{t). 
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2.1.3. Scalar higher order equations 

Scalar higher order equations can also be reduced to pseudo-differential first 
order systems. We consider 


m —1 

(2.13) 

k—0 \ct\<m—k 


E afe.o(t)D,^D> = 0 


endowed with suitable initial conditions 

B^u{0,-)=Uk, A: = 0,l,...,m-1, (2.14) 

and for coefficient functions ak,a G T{m — k— |q;|}. Assuming that the principal 
part is uniformly strictly hyperbolic, i.e., that the polynomial 

m —1 

^'” + E E = o (2.15) 

k—O |Q;|=m —fc 

has distinct real roots Xj{t, ^) uniform in t > 0 and ^ G and that a suitable 
replacement of (2.5) and of (2.7) holds true, we can provide an asymptotic 
construction of solutions. The large-time principal symbol involves only the 
asymptotic behaviour of the terms ak,oit)- 

2.2. General strategy 

To cover all the above examples it is sufficient to consider first order sys¬ 
tems with general Fourier multipliers as coefficients. Our strategy is as follows. 
We consider transformations of such systems to make them suitable for asymp¬ 
totic integration arguments. Such transformations are done in Fourier space / 
extended phase space and depend heavily on the interplay between time and 
spatial frequency. For large frequencies we apply a standard hyperbolic theory 
diagonalising the full symbol of the operator and constructing a WKB repre¬ 
sentation of the fundamental solution. The core of this strategy was given by 
K. Yagdjian in [29] and further developed in [16, 17, 18] and [20, 22] for the 
treatment of the large-time behaviour of uniformly strictly hyperbolic systems. 
Our approach for small frequencies extends [15]. We rewrite the system in 
Fuchs type and apply asymptotic integration techniques going back to Levinson 
or Hartmann-Wintner. 

2.3. Definitions: Zones, Symbol classes 

We decompose the phase space R+ x K." into zones 

Zpd(iV) = {(t,e) : (l + t)|e|<iV}, 

Z^y^{N) = {it,0 ■■ {l + m\>N}. > 

As outlined above, our strategy depends heavily on the zone under consider¬ 
ation. We denote by > 0 the function defined by (1 -I- t^)\£^\ = N, i.e.. 
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= 7V/I5I — 1 for 0 < 1^1 < iV and = 0 for |^| > N. It parameterises the 
boundary between the zones. 

To formulate our main assumptions precisely and to carry out the diagonal- 
isation procedures we need symbol classes. For this we follow [20] and define 
5{mi,m 2 } to be 

{ / 1 \ m2+fc 

aG C“(Zhyp(iV)) : |DfD^a(t,OI <k,c,N 

(2.17) 

i.e., the symbolic estimates should be satisfied on some hyperbolic zone Zhyp{N) 
for N chosen large enough. These symbol classes have natural embedding prop¬ 
erties, most important for us is 

5{mi,m 2 } C 5{m},m^l (2.18) 

whenever mi < m} and mi -I- m 2 > m} -I- m^, and behave well under multipli¬ 
cation and differentiation. We also define the residual class 

'H{k} = Pi 5{mi,m2}. (2.19) 

It plays a role as remainder class when performing perfect diagonalisation within 
the hyperbolic zone Zhyp(iV). Furthermore, symbols from the class 5{—1,2} 
are uniformly integrable with respect to t over Zhyp(iV). For a G 5{—1,2} the 
symbol estimate implies for j^j < N 

and for j^j > N 

( 2 . 21 ) 


2.4- Main assumptions 

We consider the Cauchy problem for a system (2.2), where the symbol A{t, ^) 
of the Fourier multiplier A{t,Dx) satishes the following assumptions: 

(Al) There exists a positively homogeneous symbol Ahom(CC) G i^lljO} 0 
i.e., we assume Ahom(CdC) = P^hom(CC) for d > 0, such that 

A - dlhom G 5{0,1} 0 (2.22) 

and that the eigenvalues j = 1,..., d of Ahom(C 0 s-re real and 

satisfy 

(2.23) 

uniformly in t > 0 and for all i ^ j. The system (2.2) is uniformly 
strictly hyperbolic if this assumption is satisfied. We refer to Ahom as the 
hyperbolic principal symbol. 
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(A2) Let M G 5{0,0} 0 be a diagonaliser of the hyperbolic principal 
symbol Ahom with M~^ G 5{0, 0} 0 (which exists by (Al) as shown 
in [22, Lem. 4.2]) and let 

Fo = diag {M-\A - Ai,om)M + (DtM-^)M) mod 5{-l, 2} 0 

(2.24) 

We refer to Fq G 5{0, 1} 0 as the hyperbolic subprincipal symbol. 
Then we assume that for constants G K and all t > s > 

K+ log + (7+ < Im [ Fo{t, 0 dr < K- log + C- (2.25) 
1 + s Js 1 + s 

holds true with suitable C± G M. 

(A3) There exists A G T{0} 0 diagonal and M G T{0} 0 invertible 
with M~^ G T{0} 0 such that 

R = AM - K - {tT)tM-^)M (2.26) 

is small in the sense that 

/■*« ~ dt 

sup / ||i?(t,e)|r-<oo (2.27) 


holds true for some constant cr > 1 and a zone constant N. We refer to 
to the pair (A, M) as the large time principal symbol and assume further 
that its diagonal entries A = diag(/ri,..., p,d) satisfy 

• in the case ct = 1 the weak dichotomy condition 

A d7" 

limsuplm / (/^^(t) — Pj{T)) - < oo 

or (2.28) 

cIt” 

> -oo; 


lim inf Im 

t —¥00 


• and in the case ct > 1 the strong dichotomy condition 

\lm{p.i{t) - p,j{t))\> 6 > 0, i^j. (2.29) 

(A4) For sufficiently large N the the estimate 

||DfDfA(t,e)|| S.a.jv (l + t)-i-'=+l“l (2.30) 

holds true within Zp^iN) and for all k and all multi-indices a. 

Remark 2.1. We remark that the constant N in assumptions (A2), (A3) and 
(A4) depends on the estimate we have in mind, it is determined by the number 
of diagonalisation steps to be carried out within the hyperbolic zone and a 
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sufficiently large N ensures the invertibility of certain multipliers. If (A2) is 
valid in some hyperbolic zone, it is also valid in smaller hyperbolic zones with 
possibly smaller difference k+ — k_. If (A3) is valid for some zone constant N 
it is automatically valid for all zone constants N. This follows from the uniform 
boundedness of ||.R(t,^)|| on every strip < (1 + 1)\^\ < N combined with 





(2.31) 


For practical applications, see e.g. in [15], assumption (A4) is satisfied for all 
choices of the zone constant N. 

Remark 2.2. Assumption (A2) with k_ = = 0 reduces to the generalised 

energy conservation property (GECL) introduced in [11] and also considered in 
[22]. It turns out that in this particular case uniform lower and upper bounds for 
the energy of solutions in terms of the initial energy follow, provided the Fourier 
support of the initial data does not contain = 0. On the other hand, assump¬ 
tion (Al) implies that (A2) is valid for some numbers k± and we introduce this 
assumption mainly to fix these important constants. 

Remark 2.3. If we choose in (A2) a different diagonaliser M G 5{0,0} 0 
of Ahom with M~^ G 5{0,0} 0 such that M~^Ahom-^ = Af“^AhomAf, 
then both are related by a smooth diagonal matrix JI with uniformly bounded 
entries and uniformly bounded inverse 

(2.32) 


In consequence Rg is changed to 

Fo(t,0=Fo(t,0+Dtlogff(t,0 (2.33) 

and estimate (2.25) holds true with the same constants k±. 

Remark 2.4. For some applications it is convenient to use the following stronger 
form of assumption (A3). If there exists a diagonalisable matrix Aqo G 
such that 

sup / ||fA(f,5) - Aooir—< oo (2.34) 

holds true with some exponent a > 1 and for sufficiently large N, then (A3) 
follows with the (constant) diagonaliser of Aoo as M and the resulting diagonal 
matrix as A. The weak dichotomy condition is automatically satisfied. In 
order to guarantee the strong dichotomy condition we have to assume that 
the eigenvalues of Aoo are simple. 

Remark 2.5. It is possible to generalise Assumption (A3) to matrices A and M 
which also depend on directions C/IC|. This would not change any of the results 
in Section 3.1, except for complicating notation. As this is not needed for our 
applications to pseudodifferential hyperbolic systems of first order, we omit this. 



3. Construction of fundamental solutions 


We concentrate mainly on the construction in Zpd{N). It follows [15] and 
uses in particular the tools provided in the appendix to that paper. For the 
hyperbolic zone Zhyp(iV) we will remain a bit sketchy and refer to [20] and [22] 
for the detailed procedure to be carried out. 

If U solves (2.2), its spatial Fourier transform U{t,^) satisfies the parameter 
dependent ordinary differential equation 

0 = A{t, 0U{t, 0, u{0,0 = UoiO, (3.1) 

and can thus be represented in terms of the fundamental solution S{t, s, ^), i.e., 
the solution to the matrix-valued problem 

DtS{t,s,0=A{t,0£{t,s,0, (3.2) 

as U{t,^) = £{t,0,^)Uo. Our aim is to provide an asymptotic construction of 
this fundamental solution and to derive estimates from that. 


3.1. Asymptotic integration for small frequencies 

The large time principal symbol A(t) = diag(^i(t),..., p,d{t)) € T{0}(8)C'^^‘^ 
given by assumption (A3) determines the large-time behaviour of solutions to 

m^U{t,f)=tA{t,OU{t,f) (3.3) 

within Zpd{N) as t —>■ oo. We assume further for the moment that ct = 1. Then 
the transformation [/0)(t^^) = M~^(t)U{t,^) yields 

= A(t)t/(°) -f R{t, (t, f) (3.4) 


and Levinsons theorem in the form [15, Thm. A.l] allows to construct a funda¬ 
mental system of solutions of the form 

+o(l))exp ^i^ /Xj(t)-^^ , j = l,...,d, (3.5) 

with M{t) = {vi{t)\v2{t)\ ■ ■ ’Ivdit)) G T{0} 0 having the columns Vj{t). 

As these vectors and the yij are uniformly bounded, the exponential term shows 
polynomial behaviour. This behaviour can be estimated in terms of the number 


/i = min liminf 


logt 


(3.6) 


Up to a small loss in the order this number determines the estimates also in the 
case CT > 1. 
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Theorem 3.1. Assume (AS) and let fj, be defined by (3.6). Then uniform in 
(s,^), (t,£,) G Zpd{N) with s <t the estimate 




e,N 



(3.7) 


holds true with e = 0 in the case a = 1 and for arbitrary e> 0 in the case a > 1. 

Proof. The case a = 1 just follows [15, Thm. A.l] and the two remarks after it. 
The transformed unknown = M~^U solves (3.4) with diagonal A(t) sat¬ 
isfying the (weak) dichotomy condition (2.28) and remainder R{t,() satisfying 
the integrability condition 


- dt 

sup / ||i?(t,Oll— < oo. 

?:t£>l Jl 


t 


Therefore, we find asymptotic solutions to (3.4) of the form 
{ej + o(l)) exp ^i^ , t^o 


(3.8) 


(3.9) 


uniform in ( with (t,C G 2pd{N) and with ej the j-th basis vector of C'*. 
Transforming back gives a fundamental system Vj{t,^) of (3.3) satisfying (3.5). 
The Wronskian of this fundamental system satisfies 


yVv,,...vAt) =det {Vi{t,()\---\Vdit,^)) =exp 



(3.10) 


by Liouville theorem combined with the above asymptotics. Hence, the funda¬ 
mental matrix £{t, 1, () given as 


£:(t,i,C = (U(t,CI---lCi(t,0)(n(i,CI---lCi(i,0) ' (3.11) 


can be estimated by applying Hadamard’s inequality to Cramer’s rule 


\\£{t,l,0\\<Nt-^. (3.12) 

Combined with the scaling argument from [15, Rem. A.2] 

ll^(i,ACII<iv (1^) (3.13) 

follows. For the case a > 1 equation (3.8) holds true with the exponent a. In 
order to improve integrability we apply the Hartman-Wintner theorem in the 
form [15, Thm. A.2]. This transforms the system (3.4) into a new system 

Ai {t,0 = Ait) + diag Rit, () 
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with 


sup / (3.15) 

and iteratively yields after k steps (with k chosen large enough to guarantee 
(t/ 2* < 1) a system of Levinson form. This can then be dealt with as in the first 
case. Because the new diagonal matrix Ak{t, differs from A by terms satisfying 
the cr-integrability condition, the asymptotic behaviour changes slightly. Let for 
this r(t,^) be real-valued and satisfy the u-integrability condition 

dt 

sup / |r(t,^)|'^—< oo (3.16) 

and let /i be a real number. Then by Holder inequality and with a' dual to a 
exp^^ (-Ai + ?’(T,0)-^) 

<t-^e^p(^\rU^(^j'^y ^ (3-17) 

<t-^exp(||r|U.(logt)i/"') S 

for any e > 0. □ 

3.2. Estimating derivatives 

Assumption (A4) allows to estimate derivatives of £{t^,0,^) with respect 
to the frequency variable This will be essential to prove dispersive type 
estimates. 

Theorem 3.2. Assume (AS) and (A)-). Then the estimate 

l|D|^(ic,0,e)||<a...Jv(l+t«)-'^+^|?rl“l, \^\<N (3.18) 

holds true for all multi-indices a and with arbitrary e > 0 for ct > 1 and with 
e = 0 for (7=1. 

Proof. We consider the fundamental solution £{t,0,f) for (t,^) G Zp^{N). Its 
derivatives with respect to t are easily estimated by the equation, indeed from 
Djf = A(t,f)£ we observe in combination with Theorem 3.1 that 

||D,£:(t, 0,011 < \\Ait,m ||£:(t, 0,011 S.iv (i + o-'^+^b ( 3 . 19 ) 

Higher order t-derivatives are estimated recursively taking into account the es¬ 
timate of assumption (A4). For Oderivatives we formally differentiate the equa¬ 
tion satisfied by £ and obtain 

DtD“£:(t, 0,0 = A{t, 0D^£:(t, 0,0 + E («) 0) 0,0) 

(3.20) 
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(3.21) 


such that by Duhamel formula (and using Dj£l(0, 0,^) = 0 for |q;| > 1) 


= [ S{t,s,^)Ra{s,^)ds, 
Jo 


where denotes the sum in (3.20). By induction we prove 

l|D«^(i,0,CII <o.,e,N (l + t)-''+^+l“l. 


(3.22) 


This was already shown for |a| = 0, it suffices to give the induction step. Let 
|q;| > 1. Assuming (3.22) for all multi-indices j3 < a yields in combination with 
(A4) the estimate 


Pa(i,e)ll (l + t)-''+^+l“|-b 

Hence by (3.21) in combination with Theorem 3.1 we obtain 

— fi+e 


id; 


(l+s)-A«+^+l«|-lds 

<c,e,N (l+t)-^+^ Al+s)'“'-^ds 
JO 


(3.23) 


(3.24) 


and thus the (3.22) holds true for a. Finally, mixed derivatives are easily esti¬ 
mated from differentiating (3.21). This yields 

||D,^D^“£:(t,0,CII <k,c.,e,N (l+t)-^+^-'=+l“l, lei < N, (3.25) 

and combined with the estimate 

|D^(l+tc)|<o|Cr'-'“', \i\<N, (3.26) 

the desired statement follows. □ 


3.3. The diagonalisation scheme for high frequencies 

We recall some facts about transformations to be applied within the hy¬ 
perbolic zone ^hyp(A^) in order to make the system suitable for asymptotic 
integration and WKB analysis. They are merely standard and follow [20], [22]. 
They are entirely based upon assumption (Al). First, we recall [22, Lem. 4.2]. 

Proposition 3.3. 1. The eigenvalues of the hyperbolic principal part 

Ahom{t,0 satisfy 

AjG5{1,0} (3.27) 

together with 

(A-A,)-i g5{-1,0}, ty^j. (3.28) 

2. There exists a matrix-valued .symbol M G S{0, 0} (S> homogeneous 

in the sense that M(t,pf) = M{t,f,) for p > 1, with uniformly bounded 
inverse M~^ G 5{0, 0} 0 .such that 

M~^{t, CAhom(t, OAIit, C = diag(Ai(t, f),... Ad(t, f)) (3.29) 

within 2hyp(iV). 
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This proposition allows for a first transformation of the system. If the vector¬ 
valued function U solves (2.2), the function solves 

+ Ro{t,0)v^°Ht,0 (3.30) 

with diagonal main part 

V = diag (Ai,..., Ad) G 5{1,0} (g) (3.31) 

and remainder term 

Ro = M-\A - Ahom)M + (DtM-i)M G 5{0,1} 0 (3.32) 

Note that assumption (A2) is related to the diagonal part of that symbol, 
Fq = diagi?o. Choosing the hyperbolic zone small enough allows to perform 
further transformations to improve the remainder. The first step is as follows, 
cf. Section 4.2 from [22]. 

Proposition 3.4. There exists a symbol Ni G 5{0,0} (g) satisfying G 
5{0,0} (g) for a sufficiently large zone constant N such that the operator- 
identity 

(D, - P - Ao)iVi = iVi(Di - P - Ao - Ri) (3.33) 

holds true with Fq = diagi?o G 5{1,0} (g) and Ri G 5{—1,2} (g) 

Furthermore, — I G 5{—1,1} (g) 

This can be iteratively improved. The essence is given in the following 
proposition also taken from [22]. The proof is standard and we omit it here. 

Proposition 3.5. For any number k > 1 there exists a matrix-valued symbol 
Nk G 5{0, 0} (g) satisfying G 5{0, 0} (g) for a sufficiently large 
zone constant N such that the operator-identity 


(D, - P - R^)Nk = Nk{T)t-'D- Fk-i - Rk) (3.34) 

holds true with diagonal Fk-i G 5{0,1} x satisfying F^-i — Fk -2 G 

5{1 — fc, fc} (g) and Rk G S{—k, fc -|- 1} (g) Furthermore Nk — Nk-i G 

S{-k,k}(g)C<^^‘^. 

This transformation is sufficient to provide first uniform bounds on the fun¬ 
damental solution S{t,s,f) within the hyperbolic zone Zhyp{N). We briefly 
collect the key estimates for the symbols obtained after transformation. We fix 
the number k of diagonalisation steps and in consequence also the zone constant 
N and define t^ = max{A^/|^| — 1,0} for this number N. Because Fk-i — Fq is 
uniformly integrable over Z\iyp(N) we know that 


exp(i^ Fk-i{T,f)dT 


< 


N,k 


1 -|- t 
1 -b S 


-K+ 


t > S, 


(3.35) 
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and 


exp(-i^ dr) 


< 


N,k 


1 + t 
1 + S 


t>S, 


(3.36) 


hold true with the numbers k± from (A3). If fc > k_ — the estimate 




dr <.k 


1, 


l?l > N, 
1^1 < N, 


(3.37) 


holds also true. This allows to compensate the polynomial behaviour of the 
fundamental solution of the diagonal part of the transformed system and implies 
the following first theorem. 

Theorem 3.6. Assume (Al) and (A2). Then for any k > k- — k+ + 1 we find 
a zone constant N such that within Z-typ(N) the fundamental solution £{t,t^,f) 
is representable as 

S{t,t^,f) = M{t,f)Nk{t,f)£k{£t^,f)Qk{t,t^,f)N-\t^,f)M-^m,0 (3.38) 


in terms of the matrices M from Proposition 3.3, the matrices Nk from Propo¬ 
sition 3-4, the diagonal matrices 


£kit,t^,f) +-Ffe_i(r,0) drj (3.39) 

satisfying the uniform bound 

\\£k(t,t^,f)\\ <N,k (3-40) 

and matrices Qk{t, s, f) uniformly bounded and uniformly invertible within Zhyp(A^) 
satisfying the symbolic estimates 

||D^Qfc(t,t4,e)|| ler, t>t^ (3.41) 

for all multiindices |a| < k — 1 k+ — K-. 

Combining the uniform bounds of the diagonalisers in (3.38) and using the 
estimate for with |a| = 0 we obtain in particular a uniform bound for the 
fundamental solution. 

Corollary 3.7. The estimate 

\\£it,k,m<N(^^^ (3.42) 


holds true. 
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Proof of Theorem 3.6. We construct the fundamental solution in several steps. 
First, we solve the diagonal part Dj — 2? — F^-i of the transformed system. This 
can be done explicitly and yields the fundamental matrix 

Skit,s,f) =e^p(^i (T>(t,^)+F fe_i(T,0) dr^ , (3.43) 

which can be estimated together with its inverse by (3.35) and (3.36). To 
construct the fundamental solution to the transformed system Dj — 

V — Fk-i — Rk we use the ansatz £{t, s, ^) = £k{t, s, £,)Qk{t, s, ^), which yields 

T>tQk{t,s,f) = {£k{s,t, f)Rk{t, ^)£k{t, s, f)) Qk{t, s, f) = 'jik{t,s,f)Qk{t,s,^) 

(3.44) 

together with the initial condition Qk{s,s,f) = I. The solution to this system 
can be represented in terms of the Peano-Baker series 


Qfc(f,s,0 = I + ^iW 7^fe(^l,s,0 / 7^fe(^2,s,?) 


e=i 


r^i-1 


TZk{ti,s,£,)dti-■ ■ dt 2 dti. (3.45) 


Chosing k large enough such that k > K- — k.\. implies uniform integrability of 
over the hyperbolic zone and thus the uniform bound 


||Qfc(i,f?,'?)|| < exp ||7^fe(r,^{,^)|| drj 

^ ^ dr \ C' \ 

“ 1^(1 + t^)«--«+ _^T-)fc+l-K-+K+ j - J ' 

(3.46) 


follows. It further implies the convergence Qk{t,t^,f) —>• Qk{oo,t^,f) locally 
uniform with respect to ^ 0 as well as the invertibility of Qfc(oo,t^,5) based 

on the estimate for its determinant 


det Qfc(oo, tj, ^) = exp [ i / trace7?,^(r, ^) dr 


Pe^py-dJ ||7^fc(r,7^,011 drj . (3.47) 

The fundamental solution to the original system is obtained from the diagonal- 
isation procedure. Indeed, tracing back the transformations yields (3.38). It 
remains to prove the symbolic estimate for the matrix Qk(t,t^, f). This follows 
by differentiating the series representation (3.45) term by term. To establish 
bounds for derivatives of TZk{t,s,^), we first observe that for |q;| = I 


D^£kit,s,^) = £k{t,s,£,) / 5|*(T>(r,0+i^fc_i(r,0) dr. 


(3.48) 
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where the integral on the right is bounded by (t—s). Forming further derivatives 
and using Leibniz rule implies 


||D^£:fe(t,s,^)|| t>s, (3.49) 

for all multi-indices a. Derivatives with respect to s yield 

Djk(t, S, 0 = -£k{t, S, 0 {V{s, 0 + Ffc_i (s, 0) (3.50) 

and thus are estimated by multiplications by |^|. Again using Leibniz rule and 
combining it with the derivatives with respect to ^ yields 

||DfDfgfc(t,s,OII <Ar.fc,a.^ (^1^) t>s. (3.51) 

Similarly we obtain 

||DfD“ffc(s,t,e)|| (^1^) t>s. (3.52) 

and hence estimates for derivatives of TZk{t,t^,^) follow from the symbolic be¬ 
haviour of Rk{t,^) combined with the estimates (3.51) and (3.52) and read as 

(3.53) 

Hence, differentiating (3.45) term by term with respect to ^ yields in combina¬ 
tion with (3.26) 

<iv.fc,a Icr'“', t>t^, (3.54) 

as long as |q;| < k—l — K- + K+ in order to guarantee uniform integrability of the 
appearing terms in the series and thus uniformity of the estimate with respect 
to t. As there are only finitely many multi-indices involved, the constants in the 

estimate can be chosen independent of a. □ 

4. Energy and dispersive type estimates 

Energy estimates 

Estimates for the norm of solutions follow directly from Theorems 3.1 
and 3.6. Let fi be defined by (3.6) and let k+ and k- denote the constants from 
(A2). If cr = 1 then the estimate 

0,011 < (l+t)-™i-+'A*} (4.1) 

holds true and therefore any solution U{t, x) of (2.2) satisfies the norm estimate 

l|C^(i,-)l|2<(l + t)-™^"+’^>||t/o||2. (4.2) 
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If cr > 1 we distinguish two cases. Either < /i, then we obtain similarly to 
the above situation 

\\u{t,-)\\2<{i+tr^+\m2 (4.3) 

while for k+ > fj, 

\\U{t,-)h<e{l+t)->^+^Uo\\2 (4.4) 

holds true for any e > 0 . 

Remark 4.1. Note that in contrast to results based on the GECL property intro¬ 
duced in [ 11 ] where one step of diagonalisation in the hyperbolic zone is sufficient 
to derive energy estimates and hence energy estimates essentially follow under a 
assumption for the coefficients, our energy estimates need a higher amount 
of smoothness depending on k- — 

4-2. Dispersive type estimates 

For dispersive type estimates the behaviour of high frequencies is of im¬ 
portance. First, we recall some basic facts and notations. Associated to the 
eigenvalues Ai(t, ^), ..., Xd{t,i) of the hyperbolic principal symbol Ahom(t,0 
we define the auxiliary functions 

(i> C) = ^ ^ (a 0 dr (4.5) 

and families of slowness surfaces 

E« = {eGM”:r?,(t,e) = l}. (4.6) 

For fixed t these surfaces are smooth (even algebraic if A{t, is polynomial 
with respect to ^) and disjoint. They determine the dispersive properties of the 
Fourier integral operators 

uo^ (4.7) 

for given amplitudes supported within Zhyp{N) and satisfying 

|D4“a(t,0l<c.|Cr'“', t>ti. (4.8) 

Details on these estimates can be found in [22] and will be recalled below. 
4-2.1. Assuming convexity 

(i) 

For the following we assume that the surfaces S) are strictly convex for suf¬ 
ficiently large t. We recall the notion of asymptotic contact indices for surfaces 
from [22, Chapter 4.7] and define 

7 as(Ep^;t —oo) = min {7 > 2 : liminf >r(Sp\ 7 ) > 0}, (4.9) 

t—¥OC 

where 

>f:(S, 7 ) = min>!r(E, 7 ;p) (4.10) 
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and >c{'E, 7 ;p) is defined locally around points p of S as follows. By translation 
and rotation we may assume that p is the origin and E is given as graph {y, h{y)) 
of a function h : fl ^ R and with h(0) =0 and h'(0) = 0. Then 

^ |5^h(pw)|p=o|- (4.11) 

Lower bounds on >^(E, 7 ;p) are quantitative measures of contact of order 7 
between E and its tangent plane in p. They have been introduced by the author 
and M. Ruzhansky in [20] and [21] based on earlier works of M. Sugimoto [24]. 

Theorem 4.2. Assume (Al) to (A4) and that all the slowness surfaees Ep^ 
are strictly eonvex for large t and denote 7 = miuj 7 as(Ep^; f —>■ 00 ). Then for 
l<p<2, pq=p + q, e>0 and with regularity r > {n — ~ 

following estimates hold true: 

1. If pL > K+ then 

\\U{t, Ollg S.e (1 + (4.12) 

If in addition Fq is independent of f then the estimate is valid for e = 0. 

2. If P ^ K+ then 

\\Uit, Ollg S.. (1 + (4.13) 

If in addition a = 1 in (A3) and Fq is independent off then the estimate 
is valid for e = 0. 

The constants k± are determined by (A2), the constant p by (3.6) from (AS). 

Proof. We show the first estimate. The solution U to (2.2) are represented 
as U{t,f) = S{t,0,f)Uo{f). We decompose U into two parts by means of two 
smooth cut-off functions = x(-^~^(l + 0 ICI) Xhyp = 1 —Xpd, where 

X G (^“(IR) is non-negative, bounded by 1 and satisfies x(s) = 1 for |s| < 1 and 
x(s) = 0 for |s| > 2. Then 


||Xpd(<,-)%-)lli= / Xpdit,f)\\U{t,f)\\df 

<C(l+t)-^+^/ l|C>o(C)l|de 

d(l-l-t)l«l<2W 

<C(l+t)-^+^-"||C/o|U 

based on the support of Xpd and the pointwise estimate of S{t,0,f) from The¬ 
orem 3.1. In consequence we obtain 

||Xpd(t,D,)[/(t,.)|U<C(l+f)-^+^-"||[/o||i. (4.15) 
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The treatment of Xhyp(i,Da:)C^ is more subtle. Theorem 3.6 together with 
Proposition 3.3 and Proposition 3.4 implies the representation 


d 

Xhyp{t,0U{t,0 = (4.16) 

with the phase functions ^j{t, from (4.5) and matrices Bj{t, ^) arising as prod¬ 
ucts of the diagonalisers, of the exponentials exp(J)^ Fk-i{T,^) dr) 

and S{t^,0,^). They satisfy 

||D^i?,(t,OII <C„l?r'“'(log(e + t))l“l. (4.17) 

The logarithmic term appears from derivatives of exp(/^* Fo(r, ^) dr) and dis¬ 
appears if Fq is independent of 

Finally, [22, Theorem 4.7] yields the desired decay estimates 

||Xhyp(t, D,)[/(t, OIU < 0(1 + 1)-'=++^-V WUoWbi, (4.18) 

in terms of a Besov norm of the initial data with regularity r = n — . Com¬ 

bination with (4.15) and interpolation with the already shown energy estimate 
yields the desired statement. 

The second statement follows similarly with one difference. Here the decay 
in the zone 2pd(A^) is slower and thus determines the final result. □ 

The second part of the result, where the zone Zpd(iV) determines the decay 
order, can be further improved by assuming moment and decay conditions on 
the data in analogy to the corresponding statement of [15]. 

4-2.2. No convexity assumption 

(i) 

Without assuming convexity for the slowness surfaces S) the resulting dis¬ 
persive properties are essentially one-dimensional. We recall the non-convex 
asymptotic contact indices 

7o,as(Sp^t -)> oo) = min {7 > 2 : liminf > 0}. (4.19) 

t—¥CC> 

where now 

7 

x^o(S,7) = min>f:o(S,7;p), >!ro(S,7;p) = sup V |5^h(pw)lp=o|- (4.20) 

pes 1-1=1^ 

The proof of the following statement is analogous to the convex situation just 
replacing the abstract result on Fourier integral operators for convex slowness 
surfaces by the corresponding result for general non-convex surfaces. Of course 
both situations can be mixed if some of the surfaces are convex. The overall 
decay is determined by the slowness surface implying the lowest decay rate. 
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Theorem 4.3. Assume (Al) to (A4) and denote byj = miiij 7 o,as(S('^^; t —>• oo) 
the minimal non-eonvex asymptotie contaet index of the slowness surfaces 
Then for 1 < p < 2, pq = p + q, e > 0 and with regularity r > (n — — |) 

the following estimates hold true. 

1. If pL > K+ then 

+ (4.21) 

If in addition Fq is independent of ^ then the estimate is also valid for 
e = 0. 

2. If P' ^ K+ then 


\\Uit, Oil, S,. (1 + (4.22) 

If in addition a = 1 in (A3) and Fq is independent off then the estimate 
is also valid for e = 0. 

The constants k± are determined by (A2), the constant p by (3.6) from (AS). 

5. Coming back to our examples 

In this final section we turn back to our motivating examples and provide 
some particular results following from the general theory. We restrict to specific 
examples to make the advantage of our approach more transparent. 

5.1. Differential symmetric hyperbolic systems 

Let Aj G T{0} 0 be self-adjoint matrices such that A{t,f) defined 

by (2.4) has uniform in t distinct real eigenvalues for all non-zero real f and 
B G r{l}(8)C^''‘^ such that 

dt 

\\tB{t)-B^r^ <oo (5.1) 

for some Boo G and some number cr > 1. Then assumptions (Al) to (A4) 
are clearly satisfied for 

n 

T>tU = Y,^jit)'Dco,U + B{t)U, UiO,-) = Uo, (5.2) 

i=i 

and our theory is applicable. It turns out that large-time estimates for so¬ 
lutions to (5.2) depend mainly on the properties of Boo. Indeed, as shown 
in [22, Lemma 4.2] and [20, Prop. 4 and Remark 5], there exists a family 
M{t,f) ofunitary matrices diagonalisingA(t,^), satisfying M, M* G 5{0,0} and 
Imdiag((DtM*)M) = 0. Hence, modulo integrable terms we obtain in the hy¬ 
perbolic zone Zhyp(A^) the representation Fo(t, 5) = t ^ dia,gM*(t,f)BooM{t,f) 
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and therefore the diagonal entries of tFo{t, must be elements of the numerical 
range 


of the matrix Boo- Hence, using 
K_l_ = minim 

and 

K- = maxim 
we obtain the estimates 

c* 


and 


exp(i^ Fo{T,^)dT 
exp^i^ ^b('r,C)dT 


< 


N 


< 


N 


: G CM|^;|| = 1} 

(5.3) 

= min spec Im B^o 

(5.4) 

1 = max spec Im B^o 

(5.5) 

a + t\-"+ 

,1 + J ’ 

(5.6) 

(i + tY- 

(i + J ’ 

(5.7) 


This yields the corresponding estimate for the fundamental solution 

in the hyperbolic zone Zhyp(iV). In the pseudo-differential zone Zpd(iV) the 

hyperbolic principal symbol Hhom(^ 5 C) plays no role as 



||iAhom(t,5)lly <C|C| 



dt < c|^|t{ < cN 


(5.8) 


and (A3) follows directly from (5.1). Hence, we obtain within Zpd{N) for cr = 1 


ll^(i,a,e)||<^ (5.9) 

with /i = min Im spec Hoc > k+ such that the resulting estimate is an energy 
decay estimate with rate . For a > 1 and fx = k+ we lose a small amount 
of decay, otherwise the same estimate holds true. 

Dispersive estimates depend in the behaviour of the matrices Aj{t). For the 
following we assume that the limits limt_>oo (t) = Aj(oo) exist. This implies 
that the (by assumption uniformly distinct and real) eigenvalues of 

A{t,C) = converge to the eigenvalues \k{C) of A(^) = Aj(oo)^j 

and therefore 

1 /■* 

'&k{t,0 = -: Afc(T,OdT ^ Afc(^), t^oo. (5.10) 

^ Jo 

This implies that the slowness surfaces converge to the surfaces 

eW = {e G M” : XkiO = !}• (5.11) 

The convergence in (5.10) is locally uniform for all ^-derivatives as long as ^ ^ 0. 

(k) 

Therefore the surfaces E) converge in C°° and in order to apply the results of 
[22, Chapter 4.7] it suffices to consider the surfaces 
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Theorem 5.1. Assume all the surfaces are strictly convex and let 

7 = max 7(SW). (5.12) 

Then solutions to (5.2) satisfy the estimate 

\\U{t,-)\\g<p,, + (5-13) 

for 1 <p<2, pq = p + q and regularity r > {n — ~ q'^’ If ^ ^ or 

K+ < p we may choose e = 0 , otherwise the estimate holds true for arbitrary 
e > 0. 


5.2. A particular symmetric hyperbolic system 

On K" we consider the particular symmetric hyperbolic system 


BtU 



i (hit) 

1 +1 \ 



(5.14) 


with real-valued coefficients 61,62 G T{0}. Examples of such coefficients could 
be oscillating like bj{t) = sin(aj log(t) + /3j) + jj for coefficients aj,fij,jj £ R, 
j = 1, 2 and are therefore more general as the ones just considered. We assume 
that ^ 

h “ ^^('r))-^ (5.15) 

is uniformly in t > 1 bounded from below or from above. Then for this system 
assumptions (Al) to (A4) are all satisfied. Indeed, with 


~ (lei *0) + 62(6)) ^ 


and the fact that the principal part is independent of t we see that in (A2) the 
hyperbolic subprincipal symbol is independent of It is given by 


Foit,0 


i 61 (f)-t 62 ( 6 ) ^ 
1 + t 2 


(5.17) 


Furthermore, in (A3) the diagonaliser M can be chosen as identity matrix and 
the large-time principal part is determined as A(t) = i(l-|-t)“^ diag( 6 i(t), 62 ( 6 )). 
The weak dichotomy condition is equivalent to the one-sided boundedness of 
(5.15). Thus (A3) holds true with a = 1. 

The decay rate in the zone 2ip^(JV) is thus determined by the minimum of 
the numbers 




Pi = liminf 

t->oo log t 


lim inf bj{t), lim sup bj{t) 


(5.18) 


while the decay in Zhyp(Af) is determined by their arithmetic mean (pi + p 2 )j2 
and thus stronger. If we denote p = min{/xi,/X 2 } then the main result of this 
paper implies the dispersive type estimate 

\\Uit,-)\\q <p (1 (5.19) 


ioT 1 < p < 2, pq = p + q and r > — ^). 
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5.3. A wave equation with weak dissipation 

Wave equations with weak and scale invariant dissipation 


utt-A^u + b{t)ut = Q, ^(0, •) = uo, wt(0,-)='Ui (5.20) 

have been considered by the author in [26] and [27]. The approach of [27] allows 
for real-valued & G T{1} such that 

limsupt6(t) <1 (5-21) 

t—¥00 

holds true. The present paper allows to overcome this restriction and to consider 
arbitrary non-negative b G T{ 1 } and to derive dispersive type estimates for 
solutions. Of interest is the number 

f* b{T) dr 

u = liminf22—- (5.22) 

t^oo log t 

and considerations boil down to the question whether /x < 2 or /i > 2 . 

Rewriting the problem as system of first order in 17 = (J^lu, Ditt)^, we obtain 

^(i,e)=(|^| '^')+ife(t)(o G5{1,0}+5{0,1} (5.23) 

and assumptions (Al) to (A4) are satisfied with To(f)^) = ib{t)/ 21 independent 
of M = I and for <7 = 1. The weak dichotomy condition is clearly satisfied 
for 6 > 0 and the behaviour in Zpd{N) is determined from the numbers 0 and 
fj,, while in Zhyp(iV) the overall decay stems from /x/ 2 . 

We are not looking for estimates of U in terms of Uo, but rather of norms 
for Vm and dtu in terms of Sobolev norms for ug and ui. This gives for the first 
component and small frequencies a further factor or |^| on the Fourier side (cf. 
formula (3.11) for t > 1) 

foil) O(t-n)fOil) 0(|?|)\pe|/(m (^24) 

lyC>(l) 0(1-'') j 1^0(1^]) 0(1) y ly 1 J in^pdflVj (5.24j 

= (Vl(t,«)|V2(t.«)) =(Vi(1.?)A2(1.?))-i 

and hence additional decay of one order in the first component. We distinguish 
two cases. If /i < 2, then the estimate 

||VM(l,-)||q -h l]5fu(l,-)||, <g (l+t)~^ + IlMljjp.r) (5.25) 

holds true for 1 < p < 2, pq = p + q and with r = + ^). The decay is 

determined by the behaviour of large frequencies. If /i > 2, then the estimate 

l|V'w(l,-)||g -h \\dtuit,-)\\q <q (1-71)"^ 2“^r-«^(||uollp,r+l + ||ui||p,r) (5.26) 

holds true under the same conditions and the decay is determined by the be¬ 
haviour of small frequencies. 
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5.4- Klein-Gordon equations with time-dependent potential 

Next we consider Klein-Gordon equations with t-dependent mass term, 


utt - Au + 


m?{t) 

(TTtF 


u = 0, 


u(0,-)=uo, ut(0,-)=ui. 


(5.27) 


We assume that m G T{0} is strictly non-negative. The scale-invariant case 
m(t) = mo has been treated by C. Bohme and M. Reissig in [3] by means of 
special functions. Faster decaying coefficients with (t >->■ m(t)/(11)) G L^(IR+) 
appear in C. Bohme’s PhD thesis [2], while coefficients with slower decay rate 
appear in the treatment of [4]. Our approach allows us to provide results for 
perturbations of the scale-invariant case. 

To be precise, we assume that m G T{0} satisfies 


J \m(t) - mol — < oo 


for some mo > 0 and a number a >1. As example we could consider 

OO 

m(t) = mo -b aj‘il)(2~H) 
j=o 


(5.28) 


(5.29) 


for some G C“(R) with suppy C [1,2] subject to f^ipis)-^ = 1 and a 
sequence (aj) G £'^. We apply a partial Fourier transform and rewrite (5.27) as 
system of first order in the unknown 

(5.30) 


such that 


DtU(t,0 = 


( (i+tp|£b-ri VICP + 1/(1 + 


(i+t)fiCk+i 


) 


Vvi«b+i/(i+t)" 

The coefficient matrix A(t,^) belongs to 5{1,0}. This follows from 


uit,o = Ait,ou(t,0- 

(5.31) 


|$P+m"(f) ^ ^|^|2 + i/(i+i)2g5{i 0} 

VI?P + l/(l+t)2 ^ 


and 


Using that 


i/(l -b t) 


(l+me + l 


G 5{-2,3} C5{0,1}. 


(5.32) 


(5.33) 


,_ IbP -I- 

^IF + W+^-ICI, ^^=^=U==-|e|e5{-l,2}c5{0,l} (5.34) 

VKl +1/(1 + D 
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we see that the hyperbolic principal part is given by 


^hom(^5 



(5.35) 


and thus (5.33) and (5.34) imply that = 0 in (A2). Assumption (A3) is 
satisfied with constant large-time principal symbol if we use 


A 


oo — 



(5.36) 


and it is easy to check that (A4) also holds true for any number N. In particular 
we see that spec Aqo = {i/2 ± a/tuq — 1/4}. 

The characteristic roots = ±|^| are independent of t and arise from 

the wave equation as principal part. Thus we obtain with 


M = 



mo > 5 , 

0 < Too < 5 


the dispersive type estimate 


(5.37) 


\Wt{t,-)\\q + WV Ad- 1/(1 + t)'^u{t,-)\\q 

^m,p,e (1 + i)*^ ^ 2 (p <j)(||uo||p,r+l -|- ||ui|lp_p) (5.38) 


toT 1 < p < 2, pq = p + q and r > -^^(^ — ^) generalising the results from [3]. 


5.5. Models in cosmology 

Equations with time-dependent coefficients arise in a natural way when 
studying the physics of an expanding universe. This has attracted consider¬ 
able interest over the recent years, just to mention a few references we refer to 
the work of K. Yagdjian, A. Galstian and T. Kinoshita [8, 30, 9] and sketch 
the relation to our approach. Especially for the family of Friedmann-Lemaitre- 
Robertson-Walker spacetimes on with metrics of the form 


ds^ = - dt^ -I- 


1 


n 




a{t) > 0, 


(5.39) 


in an expanding / shrinking universe the covariant Klein-Gordon equation for 
scalar fields takes the form 


2/ \ A 

Utt — a (t)Au — n—^Ut + mou = 0 
a{t) 


(5.40) 


Here and later on d denotes the t-derivative of a. Particular cases contain the 
de Sitter model with a{t) = e“*, the Einstein-de Sitter model with a{t) = 
(1 -I- or the anti-de Sitter model with a{t) = e*. A partial Liouville 
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transform allows to reduce the equation to constant coefficients in the principal 
part. Indeed, let 

A{t)= [ a{e)de (5.41) 

Jo 

and denote by r = A{t) the new time variable. Then a simple calculation yields 
for the new unknown u{t,x) = u{A~^{t),x) the equation 

Urr — + b{T)Ur + m{T)u = 0, TS[0,r], (5-42) 

with variable mass and dissipation. The time T is given by T = a{9) dO and 

may be finite. We are interested in the case when T = oo and thus excluding 
the de Sitter model. Such equations were considered in [15] and we will only 
comment on the assumptions to be made on the coefficients. The dissipation 
term is given by 


1 

II 

) 

(5.43) 

and the mass term 

, N Wo 

^Y')= 1, A-X, 

a^{A 1 r)) 

(5.44) 

In order for our approach to be applicable, we require that b G T{1} and 
that m G T{2}. This can be expressed in terms of estimates for the coeffi¬ 
cient function a(t), its derivatives and the primitive A{t). Indeed, the estimate 
|&(t)| < (7(1 -I- t)~^ is equivalent to the requirement 

|dWI / ^ a{t) 
a{t) " 1 + A{t)^ 

(5.45) 

similarly TO(r) < (7(1 -I- r)“^ reduces to 


c < j 

(5.46) 

In order to obtain all estimates for derivatives we require 
two estimates that all higher order derivatives satisfy 

in addition to these 


(5.47) 


This corresponds (in analogy) to the assumptions in [17] and immediately im¬ 
plies b G T{1} and m G T{2}. Such estimates are satisfied both for the 
Einstein-de Sitter and the anti-de Sitter model. Indeed if a{t) = (1 -I- tY for 
some number then 

A{t) = £ a{e) d0 = ^ ((1 + tY+^ - 1) , t = {l + {l+ l)r)^ , (5.48) 
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and hence 


/ 1 \ ~(^+i) f 

b{T) = (2 - n)e (^1 + (1 + (^ + l)r) j ~ (2 - n)j^T-^ + 0(r"2) 

(5.49) 

as T —>■ oo. The mass term does not have enough decay for our theory, 
m(T) = mo (l + (l + (£ + ~ mo(£+ 

(5.50) 

so we need to require mo = 0 (or else, we have to consider equations with an 
effective mass term). Similarly, we obtain for the anti-de Sitter model with 
a{t) = e* and hence A{t) = e* — 1 


6(t) 


2 — n 
1 + T ’ 


m{T) 


mo 

(l + r)2’ 


(5.51) 


which both behave in the right way and satisfy the assumptions of our theory. 
For the de Sitter model the transformation yields 6(t) = (n — 2)/(l — r) and 
m(r) = mo/(l — t)^, but restricted to the interval [0,1]. 
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